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A condensate with two internal states coupled by external electromagnetic radiation, is described 
by coupled Gross Pitaevskii equations, whose eigenstates are analogous to the dressed states of 
quantum optics. We solve for these eigenstates numerically in the case of one spatial dimension, and 
explore their properties as a function of system parameters. In contrast to the quantum optical case, 
the condensate dressed states exhibit spatial behaviour which depends on the system parameters, 
and can be manipulated by changing the cw external field. 



PACS number(s) 03.75.Fi, 05.30.Jp 
I. INTRODUCTION 

Recent experimental work with multiple species Bose 
Einstein condensates (BEC) |l],|| has motivated theoret- 
ical analysis of their wavefunctions fiEj and excitations 
Ig-Q. The primary tools used to create and investigate 
these multiple condensates are external microwave and 
radiofrcquency radiation fields, and the combined system 
is described by a set of coupled Gross Pitaevskii equations 
(GPE). In quantum optics the eigenstates of the full sys- 
tem of the (single) atom plus field are called dressed states 
(see [pi) and have proved invaluable as calculational and 
interpretational tools. In this paper we extend this con- 
cept to the eigenstates of the coupled GPE, which we 
shall call condensate dressed states. 

We consider in detail a condensate with two internal 
states |1) and |2). We solve for the eigenstates of this 
system numerically in the case of one spatial dimension, 
and with plane wave uniform intensity radiation fields. 
The major new feature that occurs in condensate dressed 
states is the spatial dependence of the wave functions, 
which can prove significant even when the external field 
is a uniform plane wave. 

We explore the properties of the condensate dressed 
states as a function of system parameters, and present re- 
sults representing broad classes of the possible behaviour. 
We begin by considering the most general properties of 
the condensate dressed states including their symmetries. 
We show that in the simplest case of identical traps 
for each component and identical collisional interactions, 
both components have identical spatial behaviour. How- 
ever, with non identical traps, the two components may 
have markedly different spatial character, and we exam- 
ine the dependence of these shapes on the trap parame- 
ters, the collisional parameters, and the external field. 
We show that the condensate's spatial shapes can be 
manipulated by changing the external field strength or 
detuning. 



II. FORMULATION 

The single particle state, |\E'), of a condensate in a su- 
perposition of the two internal states |1) and 2) can be 
written 
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where <pi(r,t) is the centre of mass meanfield wavefunc- 
tion for a particle in state \i). Under the influence of 
an electromagnetic coupling the wavefunctions associated 
with each component evolve according to the coupled 
Gross Pitaevskii equations (GPE) 
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which describe either a 1-photon M or 2-photon Raman 

f- [l2[ transition. In Eq. (||) the scaling is chosen as in 
and fi and 8l are the Rabi frequency and the detun- 
ing for the transition. The trapping potentials V\(r) and 
l^(r) for atoms in internal states |1) and |2) respectively, 
may describe any static potential, but here we restrict 
our attention to the case where they are each harmonic, 
but possibly with different spring constants and trap cen- 
tres. The eigenvalue, \x, at T = can be identified as the 
chemical potential of the system. The quantities Wu, 
w 2 2 and W12 represent the strength of the two intra- and 
the inter-species interactions, and are proportional to the 
total number of atoms within the condensed system and 
the respective scattering lengths. The state |\t) is nor- 
malised to unity, so that the fractional population n^, of 
the state \i), is given by the n, = J |-0i| 2 d 3 

We look for stationary solutions to Eq 
0i(r,t) = Vi(r)e"^*, </> 2 (r,i) = M T )e~*'' 



(§) of the form 
' which gives 



rise to the time independent coupled Gross Pitaevskii 
Equations 

M Vi = -V 2 V'i + ^i(r)^i + [w n \il>i\* + wi 2 \ip2\ 2 } i>i (5) 

n 

/iV'2 = -V 2 V'2 + V2(r)V'l + [^22 1^2 1 2 +^12^1 1 2 ] i>2 ~ (6) 
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In general this equation must be solved numerically, 
and we have used an optimisation technique (which will 
be presented elsewhere p3[) to find solutions. We present 
and discuss 

solutions which represent broad classes of possible be- 
haviour, which we discuss in the next two sections. 
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III. GENERAL PROPERTIES 

Many properties of the dressed condensate turn out 
to be analogous to those of the familiar dressed states 
of quantum optics |8|], however a key difference is the 
stationary solutions of Eq. (||) are spatially dependent. 
For simplicity we shall concentrate in this paper on solu- 
tions i\)\ and ip2 of Eq. (|5|), that have no nodes. In this 
case two eigenstates can be found, which we label \^>±) 
and denote the corresponding eigenvalues as /i±, where 
/i + > /i_. Apart from constant overall phase factors, the 
mathematical form of these states is 



|* + > = Vi+(r)|l)-lMr)|2>, 
|*_> = ^_(r)|l>+^_(r)|2), 



(8) 



where the component wavefunctions ip\± and ij>2± are 
positive real functions. Scanning the detuning of the 
electromagnetic field from the far red through to the far 
blue and solving for |\l/±) and fi± at each point reveals 
that the eigenvalues observe an avoided crossing (see Fig. 
|l|(a)), very similar to that seen in quantum optics. Here 
it is associated with a resonance in which the component 
populations are near equal (typically when \6l\ < ^), as 
shown in Fig. |l] (b). On either side of this resonance, 
as \5l\ increases, the dressed states approach a single 
component configuration (i.e the dressed state is almost 
entirely in one internal state). 

The states of Eq. (|8|) exhibit certain symmetry prop- 



erties, which relate the solution |^} = 
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parameters {£1,6l} to another solution |^) = 



for field 

1p2 

for {ft, 5l}. For the case of identical trapping potentials 
(Vi = V2) and equal collisional strengths 

(tun = w\2 — W22)) the symmetries are the same as 
those of the quantum optics dressed states, namely if the 
field is related by the rotation matrix R(6) according to 
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FIG. 1. Properties of condensate dressed states in fre- 
quency scans, (a) Eigenvalues (U+ and /i_. (b) Populations ri\ 
(solid line) and m (dashed line) of the internal states |1) and 
|2). Parameters are V1 — V2— a: 2 /4, tun = 11)12 = W22 = 500 
and Q. — 2. 
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where 



R(0) = 
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sin V cos I 
The corresponding dressed eigenvector is given by 
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with eigenvalue 
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For wn — W22 i^ w\2 the symmetry is reduced to the 
single transformation 



IJJI = "02, "02_ = Ipl- 

p, = fj,-8 L , S L = -S L . 



(13) 



IV. SPATIAL CHARACTERISTICS OF DRESSED 

STATES 



A. Effect of Different Relative Trap Potentials 

For simplicity, we will consider only harmonic traps 
with Vi — x 2 /4 and V2 = k(x — x D ) 2 /4, which allows 
condensate 2} to have a different relative spring constant 
(k) and an offset centre (x ). This arises in the JILA ex- 
periment, for example, because of the different magnetic 
moment of the two components and the effect of grav- 
ity (see H|). No exact analytic solutions are possible in 
this case, but the representative behaviour is shown in 
numerical solutions presented in Figs. and 0. In Fig. 
0(a) we see that even a small difference in the relative 
spring constant (5%) can give rise to a significant differ- 
ence in the component wavef unctions, and as k is fur- 
ther increased a clear phase separation is observed (Figs. 
0(b)- (c)), whereby one species is excluded from the region 
where the other species has high density. 

Offsets between the two trapping potentials also cause 
spatial reshaping of the two components, as can be seen 
in Fig. |3|, where x is successively increased. In Fig. 
0(a), where x is approximately 1% of the condensate 
size, significant deformation of component 2 is seen, and 
when x is increased to 1 (i.e 10% of condensate size) 
phase separation occurs (Figs. 0(b)-(c)). 



Perhaps the main interest of the dressed states is that 
they provide a simple means for manipulating the spa- 
tial shape of the condensate components. However this 
is very dependent on the relative trap potentials and col- 
lisional parameters. 

In the simplest case, where the trapping potentials are 
identical (Vi = V2) and collisional coupling coefficients 
are all equal, (wij = w), no difference in spatial charac- 
teristics of the components is possible. It is easy to show 
the that eigenstates \^±) can be written in terms of the 
eigenfunction ip of the uncoupled one component GPE 

fi V> (r) = -VVo + Vi^o(r) + w\*p o \ 2 ip (r), (14) 

and take the general form 

|*±> = ci±Vo(r) |1) T c 2 ±^o(r) |2), (15) 

with eigenvalues 
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where 
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C2± = Tp-c\ ± . (17) 

These amplitudes (cj±) and eigenvalues have precisely 
the same dependence on field parameters fl and 5l as for 
the quantum optics dressed state. 



B. Rabi Frequency Dependence 

A scan of the component wavefunction shapes as fi 
is varied is shown in Fig. [|, which illustrates the con- 
trol over the condensate profiles that is afforded by sim- 
ply altering the strength of the electromagnetic coupling 
field. Two limiting regimes can be seen. At low fields 
(0 -C \$l\) the components have distinct spatial shapes, 
determined by different trap potentials and 

collisional interactions. On the other hand, in the large 
field regime (fi ^> \6l) ^ becomes the most significant 
coupling between components and suppresses the spa- 
tial differences between the component wavef unctions. In 
particular, we notice that at large Rabi frequency phase 
separation between the condensates will be suppressed. 



C. Effect of Different Relative Collisional 
Interactions 

All of the previous results have been given for the case 
of equal collisional interactions (wn = W12 = 1022)- Al- 
though this is a good approximation to the case of Rubid- 
ium, larger variations can be expected for other atomic 
species, and it may even prove possible to manipulate the 
relative scattering lengths (e.g see |lJ|). It has previously 
been shown that the extent of phase separation of binary 
condensates in the absence of electromagnetic coupling 
depends on the relative collisional interactions ||. Here 
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FIG. 2. Effect of relative spring constant on component 
wavefunctions of dressed states for state |1) (solid line) and 
state |2) (dashed line), (a) k = 0.95, (b) k = 0.85, (c) k = 0.5. 
Parameters are as in Fig. hi, except Vi — kx 2 /4 and 8l — 2. 
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FIG. 3. Effect of trap offset on component wavefunctions. 
(a) x = 0.2, (b) x — 0.5, (c) x = 1. Parameters are as in 
Fig. 0, except V 2 = (x - x ) 2 /4 and S L - 2. 
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FIG. 4. Effect of changing Rabi frequency on component 
wavefunctions. (a) component |1). (b) component |2). The 



trapping potentials are Vi — X /4, V2 = k(x 



/4 where 



k = 



0.8 and x — 
500. 



0.5, the detuning is Sl = 2 and the 
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FIG. 5. Effect of changing wu relative to W12 and W22 on 
the component wavefunctions, with state |1) (solid line) and 
state |2) (dashed line), (a) wu = 500. (b) wu = 550. Other 
parameters are as in Fig. H (c) . 



we find analogous features arising in our dressed state 
solutions. 

In Fig. ||, we illustrate the effect of increasing wu 
relative to W22 ■ This has the effect of increasing the self 
energy of component 1, thus making it less favoured, and 
resulting in the density of component one being reduced 
in comparison to component 2. Similarly, decreasing w\\ 
makes component I energetically favoured, thus causing 
an increase in its density. 

The cross coupling term W12 mediates the interaction 
between the two components. When w\2 is larger than 
w\\ and W22, a competitive interaction occurs, so that 
the larger component in any region is favoured at the ex- 
pense of the other component, thereby resulting in den- 
sity differences being enhanced. Correspondingly, when 
W12 is smaller than wn and W22, it becomes energetically 
favourable for the components to coexist, and density dif- 
ferences are reduced. The influence that w 12 exerts on 
mixing can be clearly shown by examining how the to- 
tal populations of the two internal states change as the 
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FIG. 6. Component population (ni) in frequency scans 
with different interparticle collisional interactions. Solid line 
W12 = 500, dashed line W12 = 450 and dotted line W12 = 550. 
Other parameters are Vi = V2 — ^ 2 /4, 5l — 2, Q = 2, and 
wu = 500. 



field frequency is scanned. In Fig. |6| we see that for 
the case where all collisional rates are equal (solid line) 
the total population {fix) in state |1) changes from 0.3 
to 0.7 over the range —1.0 < 6l < 1.0. However, to see 
the same change in n\ when wi% is decreased by 10% 
(dotted line), the detuning range has to be extended to 
—2.0 < Sl < 2.0. Reducing 1012 is thus seen to favour 
mixing of the components. On the other hand, when W12 
is larger than wu, mixing is unfavourable and one state 
will usually dominate the other. This effect can be seen 
with the dashed line in Fig. o, where when the frequency 
scan passes through zero, a very abrupt reversal in the 
dominant state occurs. 



structure, in which case the individual component wave 
functions may have different numbers of nodes. In this 
context, the dressed states may be useful to help predict 
the outcome of adiabatic passage, including the possibil- 
ity of crossing to a branch of different symmetry. 
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D. Conclusion 



We have given a preliminary investigation of conden- 
sate dressed states, and explored the similarities and dif- 
ferences to the familiar dressed states of quantum optics. 
The major new feature that occurs is that the two compo- 
nents of a condensate dressed state may exhibit distinctly 
different spatial shapes, and we have examined how these 
shapes depend on the various system parameters. We 
have considered only the simplest case where the compo- 
nent wavefunctions have no nodes, and the external field 
is a uniform intensity plane wave, and shown that even 
in this case the shape of the component wavefunctions 
can be manipulated by varying the field parameters. 

Recent investigations [|HJ],[llJ have shown that when 
condensates components are coupled by a spatially vary- 
ing field, an adiabatic change of detuning may be used 
to generate an excited condensate state from the ground 
state. We can generalise our dressed state analysis to 
include the possibility of a radiation field with spatial 



